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Abstract. We investigate the behaviour of the two-point correlation function in 
the context of passive scalars for non homogeneous, non isotropic forcing ensembles. 
Exact analytical computations can be carried out in the framework of the Kraichnan 
model for each anisotropic sector. We will focus our attention on the isotropic 
sector with isotropic forcing in order to obtain a description of the influence of 
purely inhomogeneous contributions. It is shown how the homogeneous solution 
is recovered at separations smaller than an intrinsic typical lengthscale induced by 
inhomogcneities, and how the different Fourier modes in the centre-of-mass variable 
recombine themselves to give a "beating" (superposition of power laws) described by 
Bessel functions. The pure power-law behaviour is restored even if the inhomogeneous 
excitation takes place at very small scales. 
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1. Introduction 

Since the fundamental work by Kolmogorov in 1941 |I],|2] one of the key points for 
many theoretical achievements in turbulence research is the (statistical) restoration of 
homogeneity and isotropy of fluctuations at small scales [2], [I]. It is however impossible 
to provide a proper and consistent description of a great variety of systems where "non 
idealized" fluctuations are still alive. The characterization of the emission of a tracer 
from point-like sources, the study of scalar concentrations along channels with non- 
homogeneous boundaries, or in those systems whose large scales are driven by strong 
shears, are remarkable examples in which strong anisotropies and inhomogeneities 

[5] must be taken into account in order to obtain a correct picture of the statistical 
properties of those systems. 

In the last decade a consistent progress in the development of a systematic analysis 
to separate isotropic fluctuations from the anisotropic ones in real turbulent flows and 
turbulent transport 0,10], [HIE has been carried out. It has been understood how to 
treat and face systems whose rotational symmetries are broken by the presence of 
an external forcing inducing anisotropic contributions. In particular, the study of a 
simplified model for passive scalar advection, known as Kraichnan model |llj.|12j. has 
provided a clear understanding of the statistical properties in all anisotropic sectors of 
the scalar fluctuations. Indeed, closed equations for the equal-time correlation functions 
can be obtained: these are linear partial differential equations whose unforced solutions 
(also called zero modes ^3j) generally exhibit anomalous scaling. This is in contrast 
with the forced solutions that possess (non-anomalous) dimensional scaling. This has 
given the insight to explain the universality in the statistical framework. Indeed, the 
anomalous properties of small scale statistics result from a decoupling between the 
zero mode scaling and dimensional scaling, and the universality of these properties 
naturally emerges because the zero-mode scaling properties are independent on the 
forcing mechanism (see Jl] for an exhaustive review). 

In the present paper we formulate the concept of the possible small-scale homogeneity 
restoration by focusing on the two-point equal-time scalar correlation function for the 
Kraichnan advection model. The advecting velocity is still homogeneous and isotropic 
but this is not for the scalar injection mechanism which is supposed here to be neither 
isotropic nor homogeneous. As we will show, the inhomogeneous forcing induces a new 
lengthscale £ q into the scalar dynamics, in terms of which quantitative conclusions on 
the persistence of small-scale inhomogeneity will be given. The aim of this investigation 
is twofold. Firstly, we want to show how the correct homogeneous limit can be restored 
going at separations (in the two-point scalar correlation function) smaller than £ q . 
Secondly, we want to give some analytical insights about the opposite physical situation 
represented by the presence of inhomogeneous fluctuations on scales of the same order 
of our separation. In the latter regime the pure power-law behaviour (homogeneous 
limit) is replaced by the "beating" in superposition of different power laws originating 
from the scalar inhomogeneities. 
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The paper is organized as follows: in section |21 we formulate the general problem of 
anisotropic and inhomogeneous correlations. This problem will be studied analytically 
for the Kraichnan ensembles. In section El the statistical description of the isotropic 
sector is provided with particular emphasis on some kinds of physically meaningful 
forcings. Conclusions follow in section 0] 



2. The two-point non homogeneous, non isotropic scalar correlation 
functions 

2.1. Basic equations 

Let us start with the equation for a passive scalar field transported by the velocity field 
v: 

d t 6(xx,t) + v( Xl ,t) ■ V6(xi,t) = kA6(xx, t) + f(x u t) . (1) 

The velocity field v is assumed incompressible, statistically homogeneous and isotropic, 
whereas the source term / is a large-scale random forcing which is not invariant under 
translations. Let us now specialize to the case of the Kraichnan ensemble [TJ], ^2] where 
the velocity field is Gaussian, white in time, zero-mean and with two-point correlation 
function 

(v a (x 1 ,t 1 )vp(x 2 ,t 2 )) = D a/3 (xi - x 2 )6(tx - t 2 ) . 
The spatial behaviour is described by 

D aP (r) = D 5 al3 - d a p(r) , 

with 

d a p(r) = D x r 



for r = \r\ smaller than the integral scale of the velocity field (L v ), above which d a p(r) 
saturates to an almost constant value whose order of magnitude is DiL\; consequently, 
since the correlation D a p(r) has to vanish for r — > oo, the relation D ~ DiLl holds. 
d is the space dimension (> 2) and £ is the scaling exponent, describing the degree of 
roughness present in the velocity field, lying in the open interval (0,2). The forcing 
is assumed to be Gaussian, white in time, zero-mean and with two-point correlation 
function (f(x 1 ,ti)f(x 2 ,t 2 )) = F(x 1 ,x 2 )5(t 1 —t 2 ). The steady-state equation for the 
two-point equal-time scalar correlation function C(xi,x 2 ) = (6(xi,t)6(x 2 ,t)) reads 

rf - w S + s l ^ (r)+ ^ (0)l S +2K ^ + ^ +F(r ' zH0 ' (2) 

where z = (x\ + x 2 )/2 is the centre of mass and r = x\ — x 2 is the relative separation. 
To obtain (j2J) we have multiplied (JTJ) by 9(x 2 ,t), averaged over the stochastic fields 
v and / (exploiting Furutsu-Novikov's rule on integration by parts [T5]). symmetrized 
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the resulting expression and eventually dropped the temporal derivative term. Fourier 
transforming in z and defining 



C = C(r,q) = J e iq - z C(r, z)dz , 
F = F{r, q) = J e iq - z F(r, z)dz , 



we obtain, in the limit K — > 0, the corresponding equation for C: 

d a p{r) 1 f^- - \{D aP {r) + D a/3 {0)) q a q p C + F{r, q) = . 
or a orp 4 

At small separations r <C L v , since D 1 ~ D L~^, d af3 {r) ~ D (r/L v )^ is negligible with 
respect to D and D a p ~ D 5 a/ 3. In the limit r L u we can thus consider the simpler 
equation 

Generally, all the anisotropic (in the separation r) components of F(r, q) may induce 
effects that are not globally invariant under rotations. It is thus better to consider the 
following decomposition [To] : 

C(r,q) = ^2d l , m {r 1 q)Y l , m {Q) 1 

Lm 



F(r,q) = J2hm(r,q)Yi, m (ty 



Lm 



with Q denoting the solid angle associated with r, and study the behaviour of the 
correlation function at separations r smaller than the forcing correlation length L, where 
Fi jm (r, q) ~ Fi <m (0, q). The equation in each anisotropic sector (/, m) for the correlation 
function Ci iTn = C7; im (r, q) reads 

r-^dy^dA - ^±izill(d - 2 + l)r' 2 C t - £-^C t + ^ = (3) 

where, because of the degeneration, we have dropped the dependence on the subscript 
m and where we have introduced a rescaled forcing in the form (f) q j = F(0, q) / (d—l)Dx. 
A new scale £ q , associated to the strength of the scalar inhomogeneities, has also been 
introduced. It measures the separation above which inhomogeneities become important 
(e.g. it can be derived through a dimensional-analysis balance between the first and the 
third term in (J2J)) and it is defined as 

2n i -1/(2-5) 

q Do 



q [2(d-l)D 1 _ 

The Fourier scale g _1 of the forcing inhomogeneities can be larger or of the same order 
of L. In the latter case we have £ q /L ~ (L/ L v )^^ 2 ^^ <^ 1 and the lengthscale l q lies 
inside the inertial range. The homogeneous case is recovered when q~ x > L so as to 
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insure £ q /L ^> 1. An important exception to the previous limitation on g _1 is provided, 
in way of example, by the emission of a tracer from a point source, in which all g's are 
excited with the same strength. This case is however quite peculiar because L would 
be infinitesimal and £ q would fall in the diffusive range, thus preventing us from taking 
the limit k —>■ 0. In what follows we will always consider k — > and L v — > oo. 



2.2. General solution 

It is easy to verify that the general solution of (jUJ), as a function of r and £ q , is 

Cj(r; £ q ) = Q ipart (r; £ g ) + r -( d +«- 2 )/ 2 (w) + £ 2 ^M] , (4) 

where w = 2(2 - _1 (^g) (2_0/2 and v x = [(d + f - 2) 2 + 4(d + f - l)Z(d - 2 + Z)/(d - 
l)] 1 / 2 /(2 — £). The constants Bi,B 2 are fixed by the boundary conditions and K,I 
represent the Bessel functions |T2| of complex argument whose behaviours for small 
arguments are K Vl {w) ~ w~ Vl and l Vl {w) ~ wUl - The particular solution Ci ]part (r; £ q ) 
of (jnj) can be found, for instance, exploiting the method of variation of constants ^H] 
which leads to 



rw poo 

K Ul (w) p»° +1 I vl (p)dp + I Ul (w) p»° +1 K Ul (p)dp 

J0 Jw 



where A = £ q 2+d ~^ /2 ((2 - £)/2) Vo and vq = (d + ( - 2)/ (2 - £)■ Regularity at r = 
imposes -Bi = since K Vl (w) ~ while the term with i?2 is regular as r — * since 

l Vl (w) ~ w Ul . 

An exact solution can also be found for the values r>L. Indeed, in this case, we have 
Fi(r,q) ~ 0, as the forcing correlation rapidly decreases for separations much greater 
than L, and an unforced equation arises: 

r-^d r r d+ t- l dA - ^ + /_\ l) Kd - 2 + 1)^ - ^ 2 ~«Q = . (5) 
The solution of (jSJ) reduces to the zero mode, 

Q(r; £ q ) = r-(^- 2 )/ 2 [B 3 K n (w) + BJ Vl (w)} , (6) 

where regularity at oo imposes B^ = as I Vl (w) ~ w~ 1 ^ 2 e w for r — > oo, while the term 
with B 3 is regular because K Vl {w) ~ w~ 1 ^ 2 e~ w for r — > oo. 

The correlation of a passive scalar field in the presence of inhomogeneous fluctuations, 
whose characteristic length is £ q , can be thus computed in the limits of small (r <C L) 
and large (r 3> L) separations with respect to the integral scalar scale L. We find 
a dependence on some unknown constants (B 2 ,B 3 ) that can be fixed upon boundary 
conditions. To provide a clear and instructive example we can assume a forcing whose 
correlation function is a step function in r, i.e. Fi(r, q) = Fi(q)Q(L — r). In this case 
we can perform an exact matching in r = L comparing the limits (r — > L~ , r — > L + ) of 
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both Ci and C[ (prime means derivative with respect to the variable r) 
the two expressions (HJ) and © (see Appendix). The final result is 

-w 



deriving from 



Ci{r-L) 



where W 



9 2- 

'q 



+K„(w) 



P 



!'(>- 



> q ,i 



P 



^0 + 1 J 



(P)dp 



for < r < L 



for L < r < oo 



(7) 



w\ T=L = 2{2-i)-\L/Q^)l\ 
Once we assume a small-scale description with respect to the integral scale of the velocity 
field (i.e. L v — > oo), the general solution thus depends on three fundamental scales r, 
L and £ q . These scales represent the separation, the forcing correlation scale and the 
characteristic length of the inhomogeneities, respectively. From the general solution (JJJ) 
it is important to note (see Appendix) that taking the limit £ q — > oo for fixed L and 
r, (J2J) reduces to the well-known solution for the homogeneous case. Indeed, the latter 
solution satisfies the equation: 



r-^d r r d+ ^ x d r Ci 



d-l 



-l(d -2 + l)r~ 2 Ci + (j) qyl = . 



Anyway, in a small scale description with respect to the forcing correlation (r -C L), the 
presence of a finite £ g can reduce the range of pure power-law behaviour because of the 
presence of the Bessel functions in our solution. This scenario is clearly opposite to the 
one considered in the homogeneous limit where a pure power-law behaviour is found. 
In order to get a deeper insight about these two different regimes, in the next section 
we concentrate ourselves on a purely isotropic situation where the forcing correlation 
function depends only on q = \q\ and the correlation of the scalar field coincides with 
its isotropic sector (1 = 0). This is the simplest, physically relevant assumption with a 
physical meaning that one can consider to obtain a clear and systematic description of 
the influence of inhomogeneous contributions. 



3. Analysis for the isotropic case 

We focus here on the isotropic sector (I = 0) for the two-point equal-time scalar 
correlation function (the notation C is used to indicate C*/=o)- The technical advantage 
now is that alternatively to the method of variation of constants, we can choose the 
particular solution of Q to be a constant, since the second term in Q vanishes for 
I = and the coefficient of the function C(r;£ q ) reduces to a constant. We can repeat 
the same arguments given for the general case to obtain 

f^ + A 2 r^- 2 ^I u (w) forO<r<L 
lr ' q) \ A 3 r-^-^ 2 K u (w) for L < r < oo , 1 ' 

where A 2 ,A 3 are known functions of L and £ q (see Appendix). As discussed in the 
previous section, the presence of the Bessel functions (the fingerpoint of the scalar 
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inhomogeneities) makes it impossible to see a clear power-law behaviour in the inertial 
range when £ q /L ~ 1. This is clearly seen by calculating local-slopes (figure 1) of the 
difference 

S(r;£ q ) = C(0;£ q )-C(r;£ q ), 

strictly related to the second-order structure function, where C(r;£ q ) is the general 
solution (JBJ). For a fixed L and £ (say, £ = 4/3, corresponding to the KOC51 scaling 
[T9*] ) we can change the ratio £ q /L and examine the local slope behaviours as a function 
of r/L: the homogeneous case (single power- law behaviour in the inertial range) is 
recovered only when £ q ^> L, while for £ q of the same order of L a coexistence of power 
laws spoils the pure scaling of the homogeneous case. 



1.25 



1.05 



0.85 



0.65 



0.45 



0.01 



.* M. . 



0.1 

r/L 



Figure 1. Local slopes of the correlation function in the inertial range. We plot 
C( r ) — ^"dkLr^" 1 as a function of r/L for different values of the lenghtscale £ q . We plot 
IqjL = 0.1 (+), lq/L = 1 (x), l q /L = 10 (*). For comparison, the homogeneous case 
is also plotted (dotted plot). 



To be more precise, let us focus on the asymptotic properties of (JHJ). First of all we can 
perform the limit r/£ q — > to obtain 

f a(L) + b 2 (L)r 2 ^ for < r < L 
C(r; £ q ) » CU r; * 9 = V (9) 

[ h{£q)r for L < r < oo , 

where a(£ 9 ), b 2 (£ q ), b 3 (£ q ) can be obtained from the expansion of Bessel functions. In 
the limit L <C £ g , a, 62 and 63 reduce to the well-known coefficients of the homogeneous 
isotropic case a, fo, ^3 (see Appendix). 

In the opposite situation [£ q <C L) the function C(r; £ q ) approximates the step function 
CQ(L — r), where C ~ 4> q £ 2 q ~~^. Thus, if (p q = F is a constant (we shall call it "forcing of 
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the first kind"), the plot of C(r; £ q ) collapses on the axis of the abscissas when £ q — > 0. 
On the contrary, if one wanted to keep C finite, a scaling q = F £ q ^ 2 ~^ oc q 2 could be 
assumed: the collapse would now take place for £ q — > oo. This simply tells us that some 
kinds of forcing are not allowed (e.g. we also rule out 9 's giving unbounded C's or a's, 
as for example <f) q ~ V q with 7 > or < — (2 — £)). The finiteness of both a and C may 
thus be guaranteed assuming e.g. <p q = Fi£ q ^ ^ + FlL~( 2 ~& ( "forcing of the second 
kind"), with F £ , F L constants. Of course there would be infinite kinds of allowed forcing, 
but we will focus on these two because of their physical relevance. In what follows we 
will consider forcings of the first kind with F = 1 and forcings of the second kind with 
Fi — F L — 1. The value chosen for £ is its Kolmogorov value £ = 4/3 and we will focus 
on the three-dimensional case (d — 3). 

1.8 | . — . — . — . — . — . — . — . — 1 

a(£ q )/a .... 

1.4 - 



1.2 - 




Figure 2. Ratios between "actual" (functions of £ q ) and "homogeneous" (limit values 
for £ g — ► 00) coefficients: dashed lines represent a(£ g )/a, solid ones 62(^9)//? • Thin 
lines are related to the forcing of the first kind and thick lines are related to the forcing 
of the second kind. 

Figure 2 represents the ratios of the additive and multiplicative coefficients given by 
a(£ q ) and b^q) to the corresponding homogeneous ones a and /3 2 , as functions of £ q /L, 
for both kinds of forcing. The "actual" values attain the "homogeneous" ones only for 
large £ q s. 

In figure 3 we show the plots of the difference C(0;£ q ) — C(r;£ q ), together with the 
respective power-law approximations, for two different values of £ q /L, 10 2 and 10~ 2 : in 
the former case the two kinds of forcing substantially give the same result (only the first 
kind is thus represented) and the agreement is perfect all over the inertial range, while 
in the latter the separation takes place for r < £ q for both kinds of forcing. One should 
also notice that by decreasing £ q , besides the slower convergence to power-law behaviour 
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Figure 3. Plots of S(r;£ q ) — C(0;£ q ) — C(r;£ q ) and of the respective power-law 
approximations (dashed lines) for £ q /L = 10 2 and forcing of the first kind (upper plot) 
and = 10~ 2 and both kind of forcings (lower plots). Thin lines are related to the 
forcing of the first kind and thick lines are related to the forcing of the second kind. 

(as remarked in figure 1), the value C(L;£ q ) tends to decrease with the "first" kind of 
forcing and to increase with the "second" kind. 
Figure 4 shows 

AC(r L ; £ q ) = C hom (r L ; Q - C(r L ; £ q ) , 

i.e. the difference between the approximated and the actual expressions of C(r;£ q ) as 
functions of £ q /L, calculated for r = tl lying in the inertial range (in this case tl has 
been chosen as the golden section of L , but similar plots exist Vr < L). The presence 
of a maximum is quite intuitive for the first kind of forcing, as both expressions vanish 
for infinitesimal £ q , but is remarkable for the second kind, which means that the "error" 
of the approximation becomes negligible not only for large but also for small £ q s. 
The previous discussion has been carried out in the pseudo- Fourier space (r,q), but the 
final results must be expressed in the physical space (r, z) and a superposition is then 
needed. It is thus useful to analyze some instructive cases of superpositions. 
Figure 5 represents the simplest case of superposition, i.e. the excitation of two g's, 
e.g. qi and q 2 , with constant or ^-dependent amplitude. In particular the upper two 
plots are the (weighted) sum of the modes £ qi /L = 10 3 and £ q2 /L = 1 and show very 
similar behaviours between each other (departure from the r 2 ~^ straight line at r's about 
one order of magnitude smaller than the smaller £ q ) ) while in the lower two £ qi is kept 
fixed but £ q2 is reduced to 1CT 3 L. In the former case since the correlation function 
collapses towards the step function in the limit £ q — > 0, the structure function does 
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AC(r L ;£ q ) 




Figure 4. Difference between approximated (power-law behaviour with "actual" 
coefficients a(£ q ) and 62(^9)) and actual expressions of C(r) for r = golden section of L. 
The thin line is related to the forcing of the first kind while the thick one to the forcing 
of the second kind. 



2 10 




r/L 



Figure 5. Plots of the superpositions of two C(0;£ q ) — C(r;£ q ) and of the respective 
power-law approximations (dashed lines) for £ qi /L — 10 3 and £ q2 /L = 1 (upper two 
plots) or = 10~ 3 (lower two plots). Thin lines are related to the forcing of the first 
kind and thick lines are related to the forcing of the second kind. 
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not "feel" the smallest £ q in the inertial range. Obviously the restoration of the correct 
power-law behaviour depends on the degree of convergence of C(r; £ q ) towards the step 
function (in this case the first kind of forcing converges more rapidly than the second 
one). More realistic cases are connected to the excitation of a finite set of discrete 
modes: in this case the correctness of the power-law approximation is guaranteed (at 
least) for r's sufficiently smaller than the minimum £ q . On the contrary, if the forcing 
has a continuum spectrum, one has to compute the continuum Fourier antitransformed 
(J d d qe 1QZ ~ j d£ q £ q ^~^^ 2 £ q ^ d ~ 2 ^ 2 ~^^ 2 ) which is well defined for the forcings we have 
considered. 

4. Conclusions 

The properties of the two-point equal-time correlation function for the Kraichnan model 
of advection have been studied in presence of anisotropies and inhomogeneities. The 
system can be described by the following three different scales: the separation (r), the 
forcing correlation length (L) and, finally, the lengthscale of the inhomogeneities (£ q ). 
The model can be treated analytically and the properties of both small scales and large 
scales can be related to the typical lengthscale £ q . This offers the possibility to analyze 
the breaking of translationally invariant properties by means of an external forcing term 
and to check if the small scale statistics can be regarded as universal in the sense that 
it does not depend on the details of the inhomogeneous contribution. This somehow 
universal property is strictly connected to the restoration of a homogeneous limit in an 
inhomogeneous situation for scales smaller that the typical inhomogeneous one. The 
homogeneous limit (r <C £ q , L <C £ q ) has been studied and it has been shown how the 
solution reproduces exactly the one that can be obtained starting from homogeneous 
equations. On the other side, the homogeneous power-law behaviour is completely 
spoiled when £ q is of the order of the separation r and it can be seen as a "beating" of 
different power laws originating from the scalar inhomogeneities. 

Summarizing, a pure power-law behaviour exists W g going at sufficiently small r's and 
this is a clear indication of the fact that the statistical description can be seen as the 
same of the homogeneous case but with a reduced range of pure scaling law behaviour. 
When we pass to the physical space, and if more inhomogeneous modes are excited, the 
restoration of an inertial range is guaranteed if the excitation takes place only at large 
scales or at large scales together with very small scales (£ q — > 0). The calculations are 
carried out with the overall exception for those £ q s falling in the diffusive range: these 
have been excluded from our analysis. 
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We start from the equation for the projection of two-point equal-time correlation 
function C = C(r;£ q ) in the anisotropic sector I: 

r-^oy+^dA - ^±JLl11 ^ _ 2 + ly-A - £~^c l + ^ = o . 

If we assume a forcing whose correlation function is a step function in r, i.e. Fi(r, q) = 
Fi(q)Q(L — r), we can perform an exact matching in r = L by comparing the limits 
of Ci(r;£ q ) and of C[{r;£ q ) (prime means derivative respect to the variable r) deriving 
from the two expressions (jH) and ©: 

lim C^Q = C l;pait (L;£ q ) + B 2 L- {d+ ^' 2 I Ul {W) 

1 — >L~ 



lim Ci(r;4,) = B 3 L~^~ 2 ^ 2 K Ul (W) 

r—tL+ 



lim C/(r; 



d + £-2 



lim Cffr 



+ fi 2 



The correlation function must be continuous in r = L, and the same is for its first 



derivative. We can thus write the complete solution (for all £ q ) as 



Q{r;£ q ) 



Ci >pait (r;£ q ) + B 2 r-^- 2 y 2 I ul (w) 
B ^ d+ ^/ 2 K Ul (w) 
where the two constants are 



for < r < L 
for L < r < oo 



(A.1) 



B2 = -M* I p U0+1 K Ul (p)dp 
>w 



and 



B 3 = A<f> q>l / p» 0+1 K Ul (p)dp. 
Jo 

Plugging the values of B 2 ,B 3 in (|A.1|) we can obtain the exact solution written in terms 
of w, W and £ q : 



Ci(r;£ q ) 



^,i£ 2 -Zw- V0 



I Ul (w) / p»° +1 K Vl (p)dp 

J w 

p^ +l I u Xp)dp 



+K Vl {w 



^i^w-^K^w) / p V0+1 I Ul (p)dp 



for < r < L 



for L < r < oo 



In the limit £ q — > oo for fixed L and r, we find the well-known solution for the 



homogeneous case. Indeed, in this limit, for r < L we have 
'V 0+1 ^>)d/- 



(2 - 0>o - H + 2)i/j 
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and, for r > L 



P »° +1 i Vl (p)d P 



2Ki 







(2 - £) 2 (z/ + vi + 2)i/ z 



„2-£ 



\v 



(2-£)*(u + v l + 2)v l 



where (f = (-*/ - v t ){2 - 0/2 and C ; + 



C,(r) 



(2-^)2(^-^ + 2)!/, 



5 ^ 



(2-e) 2 ((^o + 2) 2 -^) 



-i/ + ^)( 2 - 0/2- We finally obtain 



" 2 € for < r < L 



for L < r < oo 



I (2-0 2 (^o + ^ + 2)^ 
which is the solution that can be exactly obained from the homogeneous equation 
projected along the anisotropic sector Z: 



r-^d r r d+ ^- 1 d r C, 



-l(d - 2 + Z)r~ 2 C 2 + 



> q ,i 







d-1 

with a step-like forcing term </> 9; /. 

In the isotropic case the particular solution of the forced equation can be chosen as 
a constant, since the second term in (JHJ) vanishes for I = and the coefficient of the 
function C(r;£ q ) reduces to a constant. The complete solution can thus be written as 
(see ©) 

A 3 r-^- 2 V 2 K u (w) 
where v stands for v\=q. The two constants are 



C{r-t q 



for < r < L 
for L < r < oo , 



A, 



and 



A, 



K V {W)^-I V {W) 



cf> q £ 2 q ^L^- 2 y 2 
K v (W)-I u (W)j 



where 



I 



d + £ 2 L -{d^)/2 Iviyv) + t (2-0/2 L -(c 2+ 2 € -2)/2 / / 



g? + e — 2 



L -(^)/2^ (1 y) + ^-(2-0/2 L -( d+ 24-2)/2 K /( M/) 



Performing the limit r, L <C ^ g and exploiting the expansion of Bessel functions, we 
easily obtain 

a(£ q ) + b 2 (£ q )r 2 -t for < r < L 



C(r;£ q 



2-d-Z 



for L < r < oo 
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where the coefficients can be found after simple but lengthy algebra: 





with 



a(£ g ) 



F(0,0)L 2 -« 



a = 



(d- 1)^(2-0(^ + ^-2) 




F(0,0) 



d(d-l)D x (2-0 




03 



F(0,0)L d 



d{d-l)D 1 (d + £-2) 



that correspond to the well-known homogeneous isotropic case. 
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